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Abstract 

Harmonic slicing has in recent years become a standard way of prescribing the 
lapse function in numerical simulations of general relativity. However, as was 
first noticed by Alcubierre [Phys. Rev. D 55, 5981 (1997)], numerical solutions 
generated using this slicing condition can show pathological behaviour. In this 
paper, analytic and numerical methods are used to examine harmonic slicings 
of Kasner and Gowdy cosmological spacetimes. It is shown that in general 
the slicings are prevented from covering the whole of the spacetimes by the 
appearance of coordinate singularities. As well as limiting the maximum 
running times of numerical simulations, the coordinate singularities can lead 
to features being produced in numerically evolved solutions which must be 
distinguished from genuine physical effects. 
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I. INTRODUCTION 



In the 3+1 formulation of general relativity, the freedom in choosing a coordinate sys- 
tem through which to describe a spacetime is replaced by freedom in choosing two gauge 
quantities: the lapse function N, which controls the slicing of the spacetime into a folia- 
tion of spatial hypersurfaces, and the shift vector N^, which defines a set of reference world 
lines threading those spatial slices. Numerical simulations based on the 3+1 formulation 
are required to prescribe a method for evaluating the lapse and the shift in terms of the 
geometrical quantities (the intrinsic metric hij and the extrinsic curvature Kij) and the 
matter fields (the density p, momentum J\ and stress S'*-') that are known on each spatial 
slice. The use of a harmonic slicing condition to determine the lapse function is appealing 
for both practical and theoretical reasons, and the incorporation of such a slicing condition 
in numerical codes is becoming increasingly common. However, recent work (Alcubierre and 
Masso , Geyer and Herold ) has found evidence that the use of harmonic slicing can 
produce pathological behaviour in numerically evolved spacetimes. 

This paper addresses the question as to how suitable harmonic slicing is for work in 
numerical relativity. In section |1| the main features of the harmonic slicing condition are 
reviewed. Following an approach similar to that of Geyer and Herold the behaviour of 
'planar' harmonic slicings of the Minkowski spacetime is analysed in section |IT|; it is demon- 
strated there that such slicings always cover the whole of the spacetime. A similar statement 
does not however hold true for the Kasner spacetime. Results derived in section |V| show 
that coordinate singularities of the type found by Alcubierre |l[] appear in harmonic slicings 
of the Kasner spacetime under reasonably general conditions. Section shows that, fur- 
thermore, the results for the Kasner spacetime carry over directly to harmonic slicings of 
the more general class of Gowdy spacetimes. In section results are presented from 
numerical simulations of the Kasner spacetime which use harmonic slicing. The coordinate 
singularities predicted by the analysis of earlier sections are indeed encountered, and the 
behaviour of the numerical solutions at times just prior to the formation of coordinate sin- 
gularities is examined. Section VH concludes the discussion by considering the implications 
that these results have for the use of harmonic slicing in numerical relativity. 



II. HARMONIC SLICING IN NUMERICAL RELATIVITY 

When considering Einstein's equations as a 3+1 evolution system, a spacetime is de- 
scribed in terms of a foliation {Ej : t G M} of spatial slices, which defines, in effect, a time 
coordinate t on the spacetime. (The details of this idea are discussed by York and famil- 
iarity with that material is assumed.) The manner in which a spacetime is sliced is a very 
important issue when the 3+1 formulation is used as a basis for performing numerical sim- 
ulations. This section recalls some basic ideas on slicing conditions in numerical relativity 
and, in particular, describes how harmonic slicings of spacetimes are constructed. 

Traditionally, work in numerical relativity has been based on maximal or constant mean 
curvature (CMC) slicing conditions, initially developed for this purpose by Eardley, Smarr 
and York among others. A spatial slice has constant mean curvature if K, the trace 

of the extrinsic curvature tensor, takes a constant value on that slice. If values for K are 
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specified across a range of slices as a function of the time coordinate t, then the 3+1 evolution 
equation for the extrinsic curvature yields an elliptic equation determining the value of the 
lapse function on each slice: 

AN - N[KijK'^ + 47r(p + S)] = -K'{t) , (1) 

where A = h^^ViVj. If the function K{t) is identically zero then the spacetime is maximally 
sliced. Maximal slicing has been found to work well in numerical simulations of asymptoti- 
cally flat spacetimes, while for closed cosmologies (in which at most one maximal slice can 
exist) the more general CMC slicing condition is appropriate. Based on their behaviour in 
simple examples and their success in numerical simulations, it is believed that the maximal 
and CMC slicing conditions will in general produce foliations which cover most, if not all, of 
a spacetime being investigated, and which at the same time avoid getting too close to any 
singularities that may form in that spacetime. 

Interest in the use of harmonic slicing has arisen in recent years because of its connec- 
tion with work that has been done in reformulating Einstein's equations as an explicitly 
hyperbolic system (see the review [§] by Reula). To date, all of the known 3+1 formulations 
of general relativity as a strongly hyperbolic evolution system with only physically rele- 
vant characteristic speeds determine the lapse through some form of the harmonic slicing 
condition. 

Harmonic slicing contrasts in several basic ways with maximal and CMC slicing. The 
lapse function in a harmonically-sliced spacetime is determined either through an alge- 
braic condition or a simple evolution equation, and not through an elliptic condition as 
in equation ([l|). This 'local' specification of the lapse is a great advantage from the point 
of view of implementing a numerical code to evolve solutions to Einstein's equations since 
elliptic partial differential equations are computationally very expensive to solve. The har- 
monic slicing condition also differs from the maximal and CMC slicing conditions in that 
it specifies a relationship between two adjacent slices in a foliation rather than determining 
properties of each individual slice: while an isolated spacetime slice can be characterized as 
maximal or of constant mean curvature, there is no such thing as an individual harmonic 
slice. One consequence of this is that the choice of lapse function on the initial slice of a 
foliation is arbitrary if the harmonic slicing condition is used. 

To simplify the following discussion, all spacetime foliations considered in this paper are 
assumed to have shift vectors A^* which are identically zero. Since the shift vector controls 
only the positioning of spatial coordinates on slices and not how the slices themselves are 
arranged, this assumption does not limit the generality of the results derived here on the 
appearance of coordinate singularities. 

The harmonic slicing condition for determining the lapse A^ can be expressed in either 
algebraic form, 

Ar = g(t,x^)v/d^t7^, (2) 

or as an evolution equation, 

^ = N{\nQ),-N'K, (3) 
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where the slicing density Q{t,x'') is an arbitrary (positive) function of the fohation coordi- 
nates which does not depend on any evolved variables. If the slicing density is independent 
of the foliation time coordinate t, then the harmonic slicing is described as simple, and in this 
case a choice of the value of the lapse N on an initial slice completely determines the value of 
the slicing density Q. The alternative situation is described as generalized harmonic slicing, 
and it is clear from equation (|^) that any spacetime foliation can be constructed using gener- 
alized harmonic slicing and a particular form for the slicing density Q. If the slicing density 
has only a 'separable' dependence on the time coordinate of the form Q{t,x^) = f{t)Q{x'^), 
then the resultant foliation has the same slices as the foliation produced by using Q{x^) as 
the slicing density, but with the slices labelled by a different time coordinate. In general it 
is not obvious how a useful slicing density Q which has a non-trivial dependence on the time 
coordinate t may be chosen, and so the simple form of harmonic slicing is most often used 
in practice. 

The main question addressed in this paper is that of how suitable the harmonic slicing 
condition is for numerical work. Several authors have already considered this question from 
various different viewpoints. Bona and Masso |^ have shown that the (simple) harmonic 
slicing condition can be used to foliate several standard spacetimes, and also that 'focusing 
singularities' are avoided by the slicing condition in much the same way as they are in 
maximally-sliced spacetimes, with the slices of the foliation not reaching the singularity in a 
finite coordinate time. (It should be noted that for a harmonically-sliced spacetime a focusing 
singularity is essentially a point at which the lapse becomes zero — this is in contrast to 
the behaviour found at the 'gauge pathologies' described below.) Cook and Scheel [|l^ have 



investigated the construction of well-behaved harmonic foliations for Kerr-Newmann black 
hole spacetimes. 

The work of Alcubierre and Masso is of particular relevance to the present discussion. 
They have shown that 'gauge pathologies' (described as 'coordinate shocks' in the earlier 
paper) can occur in numerical simulations based on hyperbolic formulations of Einstein's 
equations which use harmonic slicing. (In fact, a range of gauge conditions are considered 
by the authors, with simple harmonic slicing — the only gauge choice of interest here — 
corresponding to the special case / = 1 in their formulation. None of the alternative 
gauge choices they use correspond to the generalized harmonic slicing condition.) These 
gauge pathologies manifest as a loss of continuity at points in the evolved solution with, in 
particular, large spikes appearing in the lapse. After the time at which the pathologies appear 
the numerical solution no longer converges at the expected order. Alcubierre and Masso 
explain the appearance of gauge pathologies in terms of nonlinear behaviour in the hyperbolic 
evolution equations. However this fails to adequately answer questions about how common 
the gauge pathologies are, what happens to the foliation at the points where pathologies 
appear, and what approaches can be used to prevent the pathologies from occurring. 

Simple harmonic slicings for the Schwarzschild and Oppenheimer-Snyder spacetimes have 
been investigated by Geyer and Herold @,|]. The approach they use is based on an alter- 
native representation of equations (^ and (^ in the simple harmonic case: if T is a scalar 
function on a spacetime, the level surfaces of which represent the spatial slices of a folia- 
tion {St}, then the spacetime will be harmonically sliced (in the simple sense) if 

□T = g^^'Vf.V^T = . (4) 
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By numerically integrating equation (|^) with respect to known background metrics, Geyer 
and Herold construct simple harmonic slicings which they compare to maximal slicings of 
the spacetimes. Some properties of a slicing can be determined straightforwardly from its 
time function T, and in particular the lapse N can be evaluated through the equation 

g^''i\/,T)i\/,T) = -l/N\ (5) 

where it is clear that the vector field normal to the foliation must remain timelike if the lapse 
is to have a positive real value. In fact, Geyer and Herold find that for simple harmonic 
slicings of the Oppenheimer-Snyder spacetime, foliations which are initially timelike can at 
later times become null or spacelike, with the development of the foliation thus terminating 
at what it seems appropriate to call a coordinate singularity. The lapse becomes infinite 
as these singular points are reached, and this is consistent with the behaviour found at the 
gauge pathologies of Alcubierre and Masso. 

In what follows, Geyer and Herold's approach is applied to Minkowski, Kasner and 
Gowdy spacetimes, with the intention being to gain a better understanding of the circum- 
stances under which coordinate singularities appear, and in particular to determine whether 
they are a rare or a common feature of harmonic slicings. 



III. HARMONIC SLICINGS OF THE MINKOWSKI SPACETIME 



In the present work, equations (H) and (H) are used to investigate the formation of 
coordinate singularities in (simple) harmonic slicings of cosmological models. The following 
approach is employed. The metric Qf^^ of the spacetime being investigated is assumed to 
be known with respect to a coordinate system [t,x,y,z). An alternative foliation of the 
spacetime is constructed by taking as an initial slice one of the constant time hypersurfaces 
of the background coordinate system: the foliation time coordinate T{t, x, y, z) is given an 
initial value 



T{to,x,y,z) = to, (6) 

for some value to of the coordinate t. (It is assumed here that the hypersurface t = t^ is 
spacelike.) The lapse function of the foliation can be specified arbitrarily on the initial slice 
and this determines via equation (|^) the value of the first derivative of T away from that 
slice: 

Tt(to,x,y, z) = . (7) 

No{x,y,z)^y-g^^{to,x,y,z) 

Equations (^ and (0) provide initial data for the wave equation (Q), and the solution T 
describes a new simple harmonic slicing of the spacetime. 

As an example of how equation (^ can be used to find coordinate singularities, consider 
harmonic foliations of the Minkowski spacetime, written in standard coordinates as 

g^^dx^dx'^ = —dt^ + dx^ + dy"^ + dz^ . 
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Suppose that a 'planar' slicing of the spacetime is constructed such that the time function T 
depends only on the Minkowski coordinates t and x. Then equation takes the form of 
the one-dimensional wave equation 

T — T 

,tt -'- ,xx ; 

which has the general solution 

T{t,x) = fit + x) + h{t-x), 

for arbitrary functions f{u) and h{u). The lapse associated with this time function can be 
found from equation (^: 

= Tj -T/ 

= -Af'{t + x)h\t - x) , 

and the lapse will be well behaved as long as the function f'{t + x)h'{t — x) is positive. 

For most problems the next step in the analysis would be to use equations (^) and (|^ 
to specify a value for the lapse on an initial slice of the foliation. However the present case 
is sufficiently simple that the appearance of coordinate singularities can be studied without 
needing to specify an initial slice. If the lapse becomes infinite at a point in the foliation 
then the function f'{t + x)h'{t — x) must be zero there, and it is clear that the function must 
then be zero at all points along a line t + x = constant ot t — x = constant. Consequently, 
any spacelike slice of that foliation must include a point in it at which the lapse is infinite. It 
follows that for simple harmonic slicings of the 'planar' Minkowski spacetime no coordinate 
singularities will be present, and the foliation will cover the whole of the spacetime, provided 
that the initial slice of the foliation is everywhere spacelike. (This result is consistent with 
the numerical simulations of Minkowski spacetime reported on by Alcubierre No gauge 
pathologies are discovered for 'planar' harmonic slicings — the / = 1 case — of fiat spacetime, 
although they are found in the spherically symmetric case.) 

In the following sections a similar analysis is performed for more complicated spacetimes, 
and coordinate singularities are found to be much more common than the above result for 
the Minkowski spacetime would suggest. 

IV. HARMONIC SLICINGS OF THE KASNER SPACETIME 

In this section, equation is used to construct simple harmonic foliations of the Kasner 
spacetime, with the main point of interest being the question of whether or not coordinate 
singularities (gauge pathologies in the terminology of Alcubierre and Masso 0) form in the 
slicings. Later, in section these results are compared to numerical simulations of the 
spacetime. 

The present work uses the axisymmetric Kasner model described by the metric 

ds^ = t-^/^-dt" + dx^) + t {dy^ + dz^) , (8) 



6 



for which a cosmological singularity is present at time t = 0. For convenience in the following 
analysis (and to strengthen the connection with the Gowdy cosmologies examined in sec- 
tion |V]) a three-torus topology is imposed on the Kasner spacetime: periodicity is assumed 
in each of the coordinates x, y and z over the range [0, 27r). 

As in the simple example for the Minkowski spacetime given in the previous section, 
here 'planar' slicings of the Kasner spacetime are considered for which the time function T 
depends only on the coordinates t and x. A straightforward calculation then shows that 
for the metric @, equation for the time function of a harmonically-sliced foliation is 
equivalent to 

Ttt - T.X + = . (9) 

This equation can be solved by decomposing T spatially as a sum of Fourier modes (recalling 
that, since the model is spatially closed, the function T must be periodic in the coordinate x) 
and the general solution is found to be 

T = aolnt + bo (10) 

oo 

+ ^[Z+„(t) cos(nx) + Z_rlt) sin(nx)] , 

n=l 

with 

Z±„(t) = a±rJoi\n\t) + b±nYo{\n\t) for n 7^ , 
where a±„ and 6±„ are constants, and and Y^, are Bessel functions of the first and second 



kinds 1 1 



Equation (§) for the lapse function of the foliation can be written in the present case 

as 

Vt (T/ - T/) = , (11) 

and the requirement that the foliation be well behaved can then be expressed as 

Tt>\T,\, (12) 

where it is assumed that the orientation of the time function T of the foliation is the same as 
that of the Kasner time t, and thus that Tt is positive. If the time function T of equation (p!0D 
fails to satisfy the condition (|T2|) at a point, then the harmonic slicing must have a coordinate 
singularity there. 

The constants a±„ and 6±„ in equation (|TU|) can be determined by specifying a value A'o(a;) 
for the lapse on an initial hypersurface t = to as in equations (^ and (^. In fact, rather 
than choosing the initial value for the lapse directly, it is more convenient to choose an initial 
value I{x) for the time derivative of T, with this being specified as a Fourier sum, 

Tt{to,x) = I{x) (13) 
00 

= ko + cos(nx) -|- k^n sin(na;)] , 

n=l 
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where the constants k±n must be such that I{x) is everywhere positive. The value of the 
lapse on the initial slice can then be determined through equation (^: 

No(x) 



The particular time function T of equation ([lOD which fits the initial data of equations 
and (|I3|) is found to be 



oo 

irtr 



T = to + kotoln{t/to) ^ ^[Mnt)Yo{nto) - Yo{nt)Jo{nto)][k+nCos{nx) + A;_„sin(nx)] . 

(14) 



2 

n=l 



The simplest harmonic foliations of the Kasner metric (H) are the homogeneous ones for 
which the time function T is independent of the coordinate x. Then, to within a linear 
rescaling, the function T equals Int, and so the cosmological singularity in the spacetime 
occurs in the harmonic foliation only in the limit as T tends to negative infinity. Clearly 
the foliation covers the entire Kasner spacetime and contains no coordinate singularities. 

The behaviour of homogeneous harmonic slicings is, however, not representative of the 
general case. As an example, figure |l] shows the foliation which results from setting the lapse 
on the initial slice to be 

Noix) = I . for to = 1 , (15) 

1 + ^ sm X 

which corresponds to the choice of non-zero coefficients ko = 1 and k^i = 1/2 in equa- 
tion (1^). When equation ([Tl|) is evaluated to determine the lapse on this foliation, it is 
found that in some spacetime regions to the future of the initial slice the quantity is 
non-positive; figure |^, which plots is filled in where this happens. Coordinate sin- 

gularities must appear in any slices of the foliation which intersect these regions, and the 
development of the foliation cannot proceed beyond a time T = Tsing — 2.45 with the lapse N 
tending to infinity at points as the limiting slice is approached. In contrast, if the foliation is 
extended backwards in time from the initial slice towards the cosmological singularity, then 
the lapse appears to be well behaved on all of the slices. 

It can in fact be shown that, in general, harmonic foliations of the Kasner spacetime al- 
ways develop coordinate singularities at sufficiently late (future) times. If the time derivative 
of the function in equation (|T^) is evaluated at a fixed spatial coordinate x = Xq, then 



Ttit,Xo) = -j^ + -Y^nCnixo)[Ji{nt)Yo{nto) - Yi{nt)Jo{nto)] , 

n=l 

where it is assumed that the value Xq has been chosen such that at least one of the values c„ 
is non-zero. For sufficiently large values of the coordinate t, the approximations (9.2.5) 
and (9.2.6) of reference [|ll| can be applied to the Bessel functions in this equation. The 
result is that 
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m / \ ^0^0 / ^ \ — ^ I — \-ir / \ ( '^\T/ \ / "^M ^0^0 ^('^) 

Tt{t,Xo) ~ — — hW— to2^ V'^c„[yo('^^o)cos(nt-|7r)-Jo(nto)sm(nt-|7r)J = ^ — \ — —, 

n=l ^ 

where 0(t) is a periodic function which takes both positive and negative values. It follows 
that, for some sufficiently large value of the coordinate t, the value of must be negative 
at a point. Since this violates the condition (|12D , it must be the case that a coordinate 
singularity is present in the foliation. 

The above argument shows that coordinate singularities must in general appear in har- 
monic slicings of an expanding Kasner cosmology but gives no indication of how much of the 
spacetime a slicing will cover before it becomes pathological. To investigate this, consider 
the simple case of a time function T from equation (|14D which, like the foliation produced 
by the initial data (]T5|), contains only one non-trivial mode of amplitude k: 



T = to + htoHt/to) (16) 
- f^^[Mnt)^o{nto) - Yo{nt)Jo{nto)] sin(nx) , 

where < k < k^. For a fixed value of the Kasner time t the foliation will be well behaved 
(in that condition (|T^) is satisfied) if and only if 

° ° > max{\A{t) cos{nx)\ — B{t) sin(na;)} , 



t 



where 



A{t) = Mnt)Yo{nto) - Y^{nt)J^{nt^) 
B{t) = Ji{nt)Yointo) - Yi{nt)Jointo] 



and this condition is equivalent to 

Fit)^mtr+Bitr]<(^^J . (i?) 

It is straightforward to show that F'{t) > 0. Furthermore, if the time t is assumed to be large 



enough for approximations to be applied to the Bessel functions in A(t) and B{t), then 
F{t) ~ tn(t) where Il{t) is a positive, periodic function. An estimate of the time t = tsmg 
at which a coordinate singularity develops can then be seen to obey 

tsingoc {ko/kY , (18) 

and so, for a harmonic slicing which is perturbed from homogeneity by a single mode of 
amplitude k, the amount of the spacetime covered by the slicing becomes infinite as k tends 
to zero. 

Considering now the behaviour of the slicing in the region of spacetime between the initial 
slice and the cosmological singularity, it follows from equation (0) that, for the single mode 
case (p!6|), coordinate singularities never appear. To see this, note that the function F{t) 
must satisfy the inequality at time t = to since the lapse is required to be well behaved 
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on the initial slice, and since it is an increasing function of time it must then also satisfy 
the inequality at all earlier times. In the limit as t tends to zero (while x is constant), the 
Bessel function behaviour is known from equations (9.1.10-13) of reference |Tl|, and the 
time function T must have the form 

T ~ to[^o + kJ^ijito) cos{nx)] Int 

+ (terms that tend to a constant) . (19) 

Since the factor multiplying Int is always positive, the value of T must tend to negative 
infinity as the cosmological singularity is approached. Thus a single mode harmonic slicing 
will always cover the whole of the Kasner spacetime to the past of the initial slice. 

Returning to the general case in which an arbitrary number of modes are present in the 
time function T, the above argument can be used to derive a simple condition on the initial 



data (|T3D which ensures that a slicing is well behaved in the region of spacetime between 



the initial slice and the cosmological singularity. Suppose that equation (|l^ is split up as 

oo 

T = to + 5^ (T(+") + T(-")) , 



n=l 



where 



Ti+n) ^ k^+-)t^in(^t/to) - Kn^[Mnt)Yo{nto) - Yo{nt)Jo{nto)] cos(nx) , 
Ti-n) ^ kir'^koHt/to) - k^n'^[Mnt)Yo{nto) - Yo{nt)Mnto)] sm{nx) , 



and the values k^^^ are arbitrary subject to the condition 



oo 



n=l 



The condition ([121) that the foliation be free of coordinate singularities will then certainly 
be satisfied if 

T(+"),t > |T(+"),^| and T^-'^lt > \T'^-''\^\ Vn , 

and, as the analysis of the single mode case shows, this will be true for < t < to if it is 
true for t = to. The time derivative of each function T*^^"-* must therefore be positive on the 
initial hypersurface t = to, and for this to be the case the values k^^""^ must be chosen such 
that 

fco'^"'* > \k^n\ and k^ > \k-n\ Vn. 
This can always be done if 

oo 

ko > + \k-n\) , (20) 

n=l 
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and so any initial value for the lapse which satisfies this condition must produce a folia- 
tion which is free of coordinate singularities to the past of the initial slice. In addition to 
this, if the behaviour of the time function T is considered as t tends to zero, then it is 
straightforward to show (in analogy with equation (|19D for the single mode case) that con- 
dition (^0[) ensures that T tends to negative infinity for all values of x as the cosmological 
singularity is approached, and hence that the foliation covers all of the spacetime up to the 
singularity. The condition (^) on the initial lapse is sufficient but not necessary for the 
harmonic foliation to be well behaved to the past of the initial slice. In fact, experimenta- 
tion with different choices of initial data ( [TB| ) suggests that, even when condition ( PO] ) is not 
satisfied, harmonic foliations may never develop coordinate singularities as the cosmological 
singularity is approached, although no proof (or disproof) of this conjecture has yet been 
constructed. 



V. HARMONIC SLICINGS OF INHOMOGENEOUS COSMOLOGIES 

To summarize the results of the previous section, foliations of the Kasner spacetime (^ 
based on the simple harmonic slicing condition (^ have been investigated under the assump- 
tions that one slice of the foliation coincides with a t = constant hypersurface of the original 
metric, and that the foliation is independent of two of the standard spatial coordinates. It is 
found that when the Kasner cosmology is expanding, the slices of the foliation must eventu- 
ally develop coordinate singularities (except when the lapse is initially constant), with this 
happening at late times for foliations which are initially close to homogeneous. In contrast, 
when the Kasner cosmology is collapsing, the results suggest that coordinate singularities 
never develop in the foliation and that the slices extend all the way to the cosmological 
singularity. (This is certainly true if the lapse on the initial slice is chosen such that the 
coefficients of equation (|T^) satisfy the condition (^), and it may in fact be the case for all 
reasonable choices of initial lapse.) 

The behavioural differences found in harmonic foliations according to whether the Kas- 
ner spacetime is expanding or collapsing can be explained by considering the relationship 
between the lapse function and the mean curvature of a slice. In any foliation, the 
mean curvature provides a measure of the local convergence of world lines running normal 
to the spatial slices: 

K = -VX = -£n{ln v/diU~) , (21) 

where n'^ is the vector normal to the slices, £„ is the Lie derivative along that vector, 
and ^^/dethij is the spatial volume element. The value of K is positive when the foliation 
world lines are locally converging, and negative when they are expanding. The harmonic 
slicing condition can be written in terms of the mean curvature of the foliation using equa- 
tion (i: 

d{l/N) 

where the harmonic slicing is assumed to be simple with zero shift. It follows that when 
the mean curvature K is positive (as it may be expected to be for slices of a collapsing 
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cosmology) the value of will increase, and so the lapse will approach (but usually not 
reach) zero. (The ability of the harmonic slicing condition to avoid 'focusing singularities' at 
which the lapse vanishes is discussed by Bona and Masso [Q.) Conversely, when the mean 
curvature is negative (as it typically will be in an expanding spacetime) the value of 1/iV 
will decrease, and if it reaches zero then a coordinate singularity of the type investigated in 
this paper will develop in the foliation. 

The obvious extension of the above analysis is to the consideration of harmonic slicings 
of spacetimes more general than the Kasner model. In fact it turns out to be straightforward 
to extend the results of section |V| to a class of inhomogeneous cosmological models |]12[| : 
the unpolarized Gowdy spacetimes on (the three-torus). These spacetimes are vacuum 
models of spatially closed cosmologies with planar symmetry, and can be represented by the 
metric 

ds' = e''/H-'/^{-dt^ + de') (22) 
+ e^t[da^ + 2Q da d6 + (Q^ + e'^^) d6^' 



The spatial coordinates 6, a and 6 range from to 27r. A cosmological singularity occurs at 
time t = and the model expands forever through positive values of t. The functions A, P 
and Q depend on the coordinates 6 and r only and are required by the topology of the model 
to be periodic in 6. The vacuum Einstein equations for the metric ( p2|) reduce to evolution 
and constraint equations for the variables P, Q and A which cannot in general be solved 
exactly. Moncrief |]13| has shown that no singularities appear in the Gowdy spacetimes 



for < t < +00, and that this coordinate range describes the maximal Cauchy development 
of the models. Thus, no coordinate singularities of the type investigated in this paper can 
be present in the standard Gowdy slicing. 

If the analysis of the previous section for harmonic slicings of the axisymmetric Kasner 
spacetime (^ is applied to the Gowdy metric ( P^ ) (taking the planar symmetry of the slicing 
to coincide with the symmetry of the metric) then it is found that the time function T must 
satisfy 

Ttt - Tee + t'^Tt = , 
while the lapse N of the associated foliation is defined through 

e-^(*'^)/2Vt(T/-T/) = 1/N\ 

Comparing these expressions to equations (^ and (|TTD , it is clear that the Gowdy space- 
times admit the same harmonic foliations as the axisymmetric Kasner spacetime, and fur- 
thermore that the condition that foliations must satisfy to be free of coordinate singularities 
is the same. It thus turns out that all of the results summarized above on the behaviour of 
harmonic slicings of a Kasner cosmology carry over directly to a fairly general class of planar 
cosmological models. (In addition, it may be noted that the class of Gowdy spacetimes on 
includes as special cases the full set of Kasner spacetimes with the same topology. Thus the 



restriction of the discussion in section IV to only the axisymmetric Kasner model (H) can 



be relaxed; the same results hold for all of the Kasner spacetimes. 
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VI. NUMERICAL SIMULATIONS OF THE KASNER SPACETIME 



The results of section |rV| show that coordinate singularities are a generic feature of 
harmonic foliations of the (expanding, axisymmetric) Kasner cosmology. It follows then 
that numerical simulations of that spacetime which use the harmonic slicing condition will 
be forced to terminate after a finite number of time steps, regardless of their accuracy. The 
practical details of this are examined in the present section. 

Numerical simulations are performed based on initial data for the Kasner spacetime given 
by the metric (§) on an initial slice at time t = = 1 using the harmonic slicing condition 
and the initial value (]T5[) for the lapse. The foliation of the Kasner spacetime which results 



is the one pictured in figure |I], and the numerical simulation thus allows the prediction made 
in section ^V| — that coordinate singularities prevent the foliation from developing beyond 
a time T = Tging — 2.45 — to be verified. The slicing density of the foliation is given the 
inhomogeneous form 

f -3/4 

QiT,X,Y,Z)= \ . (23) 
1 + I sm A 

where (T, X, Y, Z) is the coordinate system of the numerical simulation, and is distinct from 
the coordinate system {t, x, y, z) of the original Kasner metric except on the initial slice 
which is labelled T = and on which (X, Y, Z) = {x, y, z). The shift vector A^' is taken to 
be identically zero, and the lack of dependence of the slicing density Q on the foliation time T 
means that the harmonic slicing is of the simple type. All of the quantities involved in the 
numerical simulation are independent of the spatial coordinates Y and Z so that the evolved 
spacetime in effect has planar symmetry. The numerical simulation is one dimensional and 
uses a high-resolution numerical scheme within an adaptive mesh refinement code to evolve a 
first-order hyperbolic formulation of the vacuum Einstein equations; a complete description 
of the code is given in reference Hl^ . 

Figure ^ displays the (T, X) coordinate system of the numerical simulation as it ap- 
pears in the (t, x) coordinates of the homogeneous background metric (^. This comparison 



between coordinate systems is made using an algorithm, described in reference |jT^, for 
tracking the positions of the observers at rest in the slices of the foliation, based on the 
values N{T, X) taken by the lapse during the course of the simulation. The observer world 
lines are plotted up until a simulation time T = 2.3. Although the inaccuracies in the 
numerical solution which are present at that time are not sufficient to stop the simulation, 
they do make it increasingly difficult to accurately track the positions of the observers. The 
T = constant surfaces reconstructed in figure |^ can be seen to be in good agreement with 
the exact results shown in figure 

The analysis of section |V| (as presented in figures |I] and |^) shows that at a time T = 
Tsing ^ 2.45 the lapse must become infinite at points of the evolved foliation. Figure ^ 
shows the actual value taken by the lapse at the earlier time T = 2.3 in the simulation, and 
it can be seen that two sharp spikes are already present in the solution. As the evolution 
progresses these spikes grow in size, and if the evolved solution were exact they would reach 
infinite heights within a short time. As they grow, the two spikes also narrow and become 
closer together. The reason for this can be seen from the paths of the coordinate observers 
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plotted in figure |^: in the region < x < vr where the coordinate singularities eventually 
form, the world lines normal to the foliation are diverging (with respect to the background 
coordinates), and so the region is resolved by a diminishing number of grid points. It 
should be noted that this divergence of observers is not itself the cause of the coordinate 
singularities; it could in principle be counteracted by an appropriate choice of shift vector for 
the foliation, but this would not prevent the lapse from becoming infinite. The narrowing 
and effective coalescence of the spikes is problematic for the simulation since features which 
are smaller than the grid spacing AX cannot be accurately resolved, and, regardless of the 
number of grid points used, eventually the numerical solution must fail to accurately model 
the behaviour of the exact solution. In fact, with the spikes being inadequately resolved, 
the evolved value of the lapse fails to become infinite (or even the computer representation 
of this) as the coordinate singularities are reached, and it is possible for the simulation to 
continue beyond the time at which the slices of the foliation cease to be spacelike in the 
exact solution. However, the numerical solution has no physical meaning past the points at 
which coordinate singularities form, and in particular, as noted by Alcubierre and Masso , 
it will no longer converge as the grid is refined. 

The spikes that develop in the lapse as the coordinate singularities are approached are at 
face value very similar to some features that are seen in numerical solutions for collapsing one- 
and two-dimensional inhomogeneous cosmologies, as studied by Hern and Stewart ||1^,T^. In 
both cases spiky features appear in the evolved variables which, as the simulation progresses, 
increase in height and decrease in width until they can no longer be resolved by the numerical 
grid used in the simulation. The question then arises as to whether the features seen in the 
inhomogeneous cosmologies have any physical relevance or whether they are, like the features 
in the Kasner simulation described above, entirely a coordinate effect. 

In reference it is concluded that the spiky features found in the inhomogeneous cos- 
mological models are not simply artifacts of the coordinate systems (or the metric variables) 
on which the simulations are based. Part of the evidence for this comes from an exami- 
nation of the profiles of the evolved variables with respect to proper distances rather than 
coordinate distances: the narrowing of the spiky features becomes even more pronounced 
when proper distances are used. It is interesting to note that the opposite result is observed 
when a similar analysis is applied to the spikes seen in the lapse in figure ^. If the metric 
component hu, which measures the proper distance in the X-direction of the simulation, 
is examined for the data plotted in figure H, it is found that sharp spikes (towards positive 
infinity) are coincident there with the spikes in the lapse, and hence that the apparently 
narrow spikes are in fact spread over large physical distances. As figure ^ demonstrates, if 
instead of being plotted against the coordinate X, the lapse is shown as a function of the 
proper distance along the X-axis, then the spikes no longer appear as notable features of the 
solution. This approach could prove to be a useful way of distinguishing genuine physical 
features in an evolved solution from effects caused by coordinate singularities in harmonic 
slicings. 
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VII. DISCUSSION 



It has been demonstrated in sections |V| and |V| of this paper that 'planar' fohations of 
Kasner and Gowdy cosmologies generated using the simple harmonic slicing condition 
will, except in special cases, always terminate at coordinate singularities without covering 
the whole of the spacetimes. When these analytic results are considered together with 
the descriptions of coordinate singularities in specific examples of 'spherically symmetric' 
harmonic foliations given by Alcubierre [|I|, and Geyer and Herold 0], a pessimistic picture 
of the nature of the harmonic slicing condition emerges: it appears likely that, except in 
special cases, foliations generated by harmonic slicing will eventually terminate at coordinate 
singularities. (Seemingly this pathological behaviour can occur whenever the foliation locally 
undergoes a protracted period of expansion). 

The coordinate singularities examined here arise as part of exact solutions for 
harmonically-sliced spacetimes, and they cannot be avoided simply by employing alternative 
numerical methods in simulations. This is clearly a drawback to the use of harmonic slic- 
ing in numerical work, and by extension to the use of hyperbolic formulations of Einstein's 
equations which typically rely on this type of slicing. It should be appreciated however that 
no known slicing condition is guaranteed to completely cover an arbitrary spacetime, and 
potential problems with coordinate singularities do not necessarily outweigh the advantages 
of using hyperbolic formulations in numerical relativity. In section ^ the behaviour of a 
numerically evolved solution is examined as it approaches a coordinate singularity, and it is 
found that by considering measures of proper distance it may be possible to distinguish the 
effects of coordinate singularities from genuine physical features of a spacetime. In general 
though the most reliable approach to identifying the presence of a coordinate singularity in 
a numerical solution is the one recognized by Alcubierre and Masso 0: convergence of the 
solution is lost when a coordinate singularity develops. 

Of course, the formation of coordinate singularities has only been discussed here for the 
simple harmonic slicing condition, and it is possible that the freedom in choosing the time 
dependence of the shcing density in generalized harmonic slicing could be put to use in 
controlling the development of the spacetime foliation such that coordinate singularities are 
avoided. Although when used in the context of a hyperbolic formulation the slicing density 
is formally required to be independent of the evolved variables, in practice there seems to 
be no problem in allowing occasional 'corrections' to be made to its value in response to 
the behaviour of the foliation; in effect this amounts to intermittently halting the simulation 
and choosing a new value for the lapse on the current slice. (An important point here is that 
changes to the slicing density should be made only at fixed time intervals, rather than after 
a fixed number of time steps, since otherwise the exact solution being sought will depend on 
the resolution of the simulation.) This 'piecewise harmonic' form of slicing has been used in 
some preliminary tests employing simple heuristics for making alterations to the value of the 
slicing density, with the resulting foliations being examined using the world line integration 
algorithm which was used to produce figure ^ (The algorithm was in fact developed for this 
purpose.) As yet however no definite conclusions regarding the effectiveness of this approach 
have been reached. 
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FIGURES 




FIG. 1. A simple harmonic foliation as seen in the background Kasner metric. The initial slice of 
the foliation coincides with the t = 1 hypersurface of the metric (js]) and is given an inhomogeneous 
initial lapse (^). The time function T is evaluated from equation (|14|) and presented as a contour 
plot with respect to the coordinates (t, x) of the original metric. The contours represent slices of 
the foliation with the dashed line marking the position of the initial slice. 
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FIG. 2. Lapse function for the harmonic foUation of figure |^. The lapse N is evaluated from 
the foliation time T through equation ( [ll| ) and presented with respect to the coordinates [t, x) of 
the background Kasner metric (^). Evenly spaced contours of the function l/N"^ are shown. The 
filled-in regions at the top of the plot show where ceases to be a positive function, and thus 

where coordinate singularities must appear in the slices of the foliation. As in figure ||, the dashed 
line marks the position of the initial slice. 
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FIG. 3. Coordinate world lines for a numerical simulation as seen in the background Kasner 
metric. A homogeneous initial data set is numerically evolved using harmonic slicing based on the 
inhomogeneous gauge function (p3|). The world lines of the simulation's grid points are tracked 
relative to the coordinate system [t, x) of the background metric (^) using the numerical values 
produced for the lapse function. Two hundred grid points are evenly spaced on the t = 1 hypersur- 
face of the background metric at the initial simulation time T = 1, and the paths they subsequently 
follow are plotted up until the end of the simulation at time T = 2.3. Periodicity of the spatial 
domain means that observers who 'go off' at x = will 'come on' at a; = 2tt. This plot is to be 
compared with figure |^. 
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lapse A' at simulation time T= 2.3 




simulation coordinate X 



FIG. 4. Lapse function for the inhomogeneously-sliced Kasner spacetime. The sohd hne shows 
the value of the lapse at the final time of the numerical simulation presented in figure ^. Co- 
ordinate singularities are responsible for the sharp spikes in the data. (The dotted line shows an 
alternative value for the lapse evaluated indirectly by combining the estimated world line positions 
of figure |^ with the analytic lapse values plotted in figure ^; the results of both approaches can be 
seen to be in good agreement.) 



lapse N at simulation time T= 2.3 




proper distance along X-axis 



5.25 



FIG. 5. The lapse from figure ^ is re-plotted against proper distance in the X-direction. The 
peaks which in the original plot appear very sharp are found to be quite smooth when viewed 
with respect to a physical measure of distance. (The proper distance is measured starting from 
the point X = 0, and an orbit of the evolved spacetime at time T = 2.3 is found to have a proper 
length of approximately 5.25.) 
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